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CHAPTER I 


INTRODUCTION 

In recent years interest has grown in the use of 
nonstationary random processes for modelling physical 
phenomena. This is due to the fact that many physical 
phenomena display nonstationary behavior. Currently these 
phenomena are generally modelled as stationary because all 
aspects of the theory of stationary processes e.g. prediction 
and spectral representation are quite well developed [1] and 
[7] . However, if these phenomena are to be better understood, 
they must be modeled as nonstationary. One aspect which 
distinguishes a nonstationary process from a stationary 
process is the power spectral density. Nonstationary 
processes have a two dimensional spectral representation while 
stationary processes have a one dimensional spectral 
representation. For nonstationary processes the theory is not 
complete, some basic questions such as interpretations of the 
spectral representation still need to be investigated. 

In this thesis, a two dimensional estimate for the power 
spectral density of a nonstationary process will be developed. 
The estimate will be applied to helicopter noise data which 
is clearly nonstationary. The acoustic pressure from the 
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isolated main rotor and isolated tail rotor is known to be 
periodically correlated (PC) and the combined noise from the 
main and tail rotors is assumed to be correlation 
autoregressive (CAR) . The results of this nonstationary 
analysis will be compared with the current method of assuming 
that the data is stationary and analyzing it as such. Another 
method of analysis is to introduce a random phase shift into 
the data as shown by Papoulis [8] to produce a time history 
which can then be accurately modeled as stationary. This 
method will also be investigated for the helicopter data in 
this thesis. 

The chapters of this thesis will be outlined as follows: 
the remainder of chapter I discusses background material 
necessary for understanding the development of this thesis. 
.Chapter II discusses a method used to determine the period of 
a PC process when the period is not known. The period of a 
PC process must be known in order to produce an accurate 
spectral representation for the process. In chapter III, the 
spectral estimate is actually developed. The bias and 
variability of the estimate are also discussed. Finally, in 
chapter IV, the current method for analyzing nonstationary 
data is compared to that of using a two dimensional spectral 
representation. In addition, the method of phase shifting the 
data is examined. Conclusions axe then made regarding the 
comparison of these methods. 
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A. Preliminaries and Background 

The purpose of this section is to review basic 
definitions and to introduce the classes of stochastic 
processes which will be investigated in this thesis. 

Definitions 

The intent here is to give a review of some fundamental 
concepts and definitions of probability theory and stochastic 
process theory. 

The probability space associated with a random experiment 
consists of three items (S,/J,P): 

1. S is the sample space containing all possible outcomes S’. 

2. j3 is the collection of all events or subsets of S. 

3. P is the probability measure defined on /3 . 

A random variable on a probability space (S,/3,P) is any 
^-measurable function X that usually maps the sample space S 
into the real line R (although the function X can take complex 
values) . Any real random variable X has a distribution 
function F x defined by 

F x (a) = P(s:eS: X(s-)<a}, for all aeR. 

The distribution function has four basic properties. 

1. F x (+°°) = 1 

2. F x ( -°°) = 0 

3. F x (a^) < F x (a 2 ) for a x <a 2 

4. F x is continuous from the right, that is lim F x (a) = 

a-+a 0 

F x ( a o ) = F x (a 0 ). 
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If it exists, there is also a density function p x 
associated with a random variable X such that 

a 

F x (a) = J*Px ( 7 ) dTf for all aeR. 

-00 

The density function p x (a)>0 for all aeR. 

The moments of a random variable X are defined by 

00 

m x k = E(X k ) = J a k p x (a) da, k=0,l,2,..., 

-00 

where k is the order of the moment. The central moments of 
X are defined as follows 

oo 

77 x k = E{(X-M) k } = J (a-M) k Px( a ) d0£ / k=0,l,2,..., 

—00 

where k is the order of the central moment. Here, the first 
order moment vx x 1= n is called the mean or expected value. And 
the second order central moment is called the variance of X. 

A stochastic process or random process is a family of 
random variables {X(A,f), AeA}, where A is the index set of 
the parameter A . Usually the index set A is either the set 
Z = { 0 , ±1, ±2 , . . . } of integers in which case the process is 
called discrete, or the set R of real numbers, in which case 
the process is called continuous. Although the developments 
made in this work are for a continuous process, the actual 
data used to produce the spectral estimate is discrete, due 
to the necessity of sampling. 

The moments of a stochastic process are defined similar 
to those of a random variable except that they now depend on 
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time. The first moment is now called the mean function, and 
the second central moment is called the variance function of 
the process. 

B. Classes of Stochastic Processes 

Let X ( t) be a stochastic process with a finite second 
moment whose mean is zero and let 

R x (t 1 ,t 2 )=E(X(t 1 )X(t 2 ) }. 

This relationship is called the autocorrelation of the 
process . 

A stochastic process X(t) is called stationary if 
R^t^t;,) depends only on Therefore, 

RxCti/tz) = R x (t!+t, t 2 +t) for all t x , t 2 eR 
where t is an arbitrary time. Stationary processes are the 
most well-known process and have been thoroughly developed, 
see for example [1]. In general any process which is not 
stationary is called nonstationary. Stationary processes are 
a subclass of the following nonstationary processes: PC, CAR, 
and harmonizable processes. All stationary processes have 
a spectral representation, a fact, which is very useful for 
obtaining information about the process. 

A process X(t) is called periodically correlated (PC) if 
there exists a time T>0 such that 

RxCt^tz) = R x (t 1 +T, t 2 +T) for all t^t^R. 

The smallest such T is called period of the process. A class 
of processes called correlation autoregressive (CAR) was 
introduced by Miamee and Hardin in [5] and [3]. CAR processes 
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are defined to be those processes for which there exist 
finitely many scalars aj such that 

R x (t 1/ t 2 ) = S Rfti+Y^ta+Yj) , for all t x ,t 2 6R 

j 

where the Y i ' s are fixed times. 

A random process X(t) is called harmonizable if the 
double Fourier transform of its autocorrelation exists. That 


IS 

00 

00 ■“ i ( <*) ^t l'”0> 2 t 2 ) 

(1) 

S x («i,w 2 ) = 1/47T 2 J dti 

J dt 2 R x (t 1 ,t 2 )e 

-00 

-00 

If 

this is the case, then 

is called the power 


spectral density of the process X(t). All of the above 
classes of processes have some members which are harmonizable 
(see Figure 1). For a stationary process, replacing R(t 2 ,t 2 ) 
by R(t 2 -t x ) and substituting t 2 -t x by r allows (1) to simplify 

to 

S x (w 1 ,w 2 ) = S x { (w 1 +w 2 ) /2 } S (u> 2 -w 1 ) 
where 6(*) is the Dirac delta function, and 

oo -io>T 

S x (w) = 1/2 jt J R *(r)e dr, 

—00 

is a one dimensional spectral density. This relationship will 
be developed in more detail in chapter III. 



CHAPTER II 


DETERMINING THE PERIOD OF A PC PROCESS 

In practice, one usually starts with a finite amount of 
discrete data. This data is generally taken with respect to 
time. According to Fourier analysis, a time history can be 
decomposed into sinusoids of different frequencies. It is 
this decomposition that we would like to use in order to 
produce an estimate of the power spectral density of a 
process . 

The Fourier transform of discrete data can be produced 
by using a computer program which can be found in the library 
of most operating systems. The discrete Fourier transform of 
a vector of sampled data {X n , n=0 , 1 , . . . ,N— 1 } , where N is the 
number of samples taken can be written as 

N-l 

X k = 2 X n e (i2 " kn/H) , 
n=0 

and {X k , k=0, 1, . . . ,N-1} goes from zero frequency to 27r/At by 
intervals of length 27r/NAt where At is the time between 
samples . 

Thus, the transformed data has the same length as the 
original data. Also it is given with respect to frequencies 
which go from the zero frequency to some maximum frequency in 
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equal intervals. It is this characteristic of the transformed 
data taking values at frequencies of equal intervals which is 
of concern when one is trying to produce the spectrum of a PC 
process . 

There is also a cutoff frequency 

— 7T/ At f 

which is the highest frequency that can be reproduced from 
data sampled at equal intervals At. This frequency is called 
the Nyquist frequency o> c . This phenomenon can be understood 
mathematically by taking a sinusoid X(t)=Acos Wl t of frequency 

Then discretizing 

= A cosujnAt 
= A cos ( w c + o 2 )nAt 
= A cos(u c nAt + w 2 nAt) 

= A cos(7rn + w 2 nAt) 

= A(cos 7 rn cosw 2 nAt — sinTrn sinw 2 nAt) 

= A cos(n7r - w 2 nAt) 

= A cos ( w c — « 2 ) nAt = A cosw a nAt 

Therefore, the frequency W] =w c +w 2 is 
indistinguishable for sampled data from the frequency 
Wa =u c -u 2 ‘ For this reason, when analyzing sampled data, X k , 
k=0, 1, . . . ,N/2 will be used corresponding to frequencies from 

zero to 7T / At . 

In this thesis, the discrete data used was taken from an 
acoustic time history produced by a helicopter fixed with 
respect to the observer. The passing blades from an isolated 
main rotor or isolated tail rotor will produce a periodic 
sound pressure time history. This time history will have the 


w 1 =w c +w 2 where w 2 <w c . 

X (nAt) 


where w a =w c -« 2 . 
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same period as that of the passing blades which produce them. 
The history produced by the main or tail rotor rotating 
alonethus represents a PC process whose correlation has the 
same period as the time history. Recall that a periodically 
correlated process is a process whose correlation is a 
periodic function. Therefore, its correlation has a period 
T and a frequency f associated with it, and it is this period 
that one must find, first. 

By definition, a function F(x) is called periodic if 
there is a positive number T such that F(x)=F(x+T) for all x. 
The smallest such T is called the period. The reciprocal of 
the period is called the fundamental frequency f 0 of the 
function, i.e. 

f 0 =_JL_ • 

T 

Frequency corresponds to the number of cycles per second and 
has units in hertz. Frequency can also be viewed in terms of 
radian frequency w which has units of radians per second and 
is related to f by 

w = 27Tf . 

The Fourier transform of discrete data taken from a 
periodic time history should only have values at frequencies 
corresponding to the fundamental frequency and integer 
multiples of the fundamental frequency [2]. However, due to 
the nature of the transform program, values are given for the 
transformed data at frequencies of equal intervals. The 
location of these frequencies will depend on the length of the 
original input data to the program. Therefore, in order to 
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avoid having values at frequencies at which there should be 
zero values, data must be input to the program of lengths 
equal to the period or some integer multiple of the period for 
a PC process [2] . 

In this chapter, a method will be examined for 

determining the period of a periodically correlated process 
from given data. The method entails finding the lines of 
support for the PC process and then to use the distance 
between these lines to determine the period of the process. 

Spectral Support 

Hurd [4] introduces a technique which is useful for 

determining the support of a nonstationary process. The 

technique consists of producing the discrete Fourier transform 
X k from sampled data of a nonstationary process. Products of 
transforms, X p X q , are then plotted in the p,q plane. Subsets 
of these products are then summed along diagonals and 

normalized with respect to the main diagonal. That is, a 
spectral coherence is produced at coordinates (p,q) where p 
and q correspond to the 27rp/NAt and 27rq/NAt frequencies 
respectively. The coherence is produced by 

M-l_ _ 

| S X p+tn X q+m | 2 
m=0 

|T(p,q,M)| 2 = 

M-l_ M-l_ 

s I x p+m 1 2 S|x q ^| 2 

m=0 m=0 

This coherence is used to determine which points over the 
array being considered have significant values. This is 
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determined by choosing a threshold value and plotting points 
(p,q) for which the coherence exceeds the threshold. If the 
process is PC, the ploted points (p,q) should produce a graph 
(Figure 5) of dots along parallel diagonal lines. The 
separation between these diagonals can be used to determine 
whether the process is PC, and to find the period. 

This technique is based on the theoretical result that 
the spectral support for a PC process with period T is on 
equidistant straight lines parallel to the main diagonal [5]. 
That is D k lines where 

D k = {(p,q): p=q+2irk} , k=0,±l, — ,±(T-1). 

T 

Taking the case where k=l, the spacing of these lines is given 
by 

p-q = 2 tt 
T 

and letting T=nAt, produces 

p-q = 27 t 

nAt . 

The difference p-q is found from the graph. Thus the number 
n of time intervals At needed to produce the period of the 
process can be calculated directly. After the period has been 
found, our next task is to estimate the power spectral 
density. 



CHAPTER III 


A POWER SPECTRAL ESTIMATION TECHNIQUE 
Ws will first consider the usual requirements which are 

A 

placed on an estimate. For a random variable A which is an 
estimate of an unknown parameter A , a usual requirement is 
that 

A 

E(A) = A, 

where the expectation is over all possible values of the 
random variable. When this is the case, the estimate is said 
to be unbiased. In addition, an estimate is chosen such that 
its variation about its mean is as small as possible. This 
uncertainty of an estimate is measured by the standard 
deviation. That is, 

A A 

Uncertainty = (E{ [ A-E(A) ] 2 } ) 1/z . 

Therefore, the usual requirements are that the estimate be 
unbiased and have the smallest possible uncertainty. 

For stationary processes, there are two well known 
techniques for spectral estimation which are the Blackman — 
Tukey and the finite Fourier transform techniques. To 
produce an estimate for the power spectral density, we have 
chosen the finite Fourier transform technique. This method 
consists of taking the discrete Fourier transform of sampled 
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data and using the transformed data to produce a spectral 
estimate. Consider a stationary process X(t) of which the 
(generalized) Fourier transform is given as 

00 — iwt 

X ( w) = _1_ J X (t) e dt 

2 7T -« 

Since X(w) is complex, its autocorrelation is given as 

CO 00 — i ( Witi — W 2 t 2 ) 

E(X( Wl )X*(w 2 )) = _1 X dt! X dt 2 RxCtz-tiJe 

47f- Z -00 -oo 

where the asterisk indicates a complex conjugate. Introducing 
the change of variables 

t = t x + t, and r=t 2 - tj 
2 

produces 

oo i ( ) t/2 

E (X (wj) X* (w 2 ) ) = X dT R K^ T ) e 

2 n -oo 

oo — i (co—o ) t 

x 1 X dt e 
2n -« 

= S x { (W!+W 2 ) /2 ) 

Therefore, 

00 

S x ( w) = X E[X(«) X* («*>') ]d« 1 
—00 

In practical situations only a single sample function of 
length T of a random process X(t) is available. Based on the 
above relationship just obtained, a class of power spectral 
estimates 

S x (co) = W.|X f («) | 2 


is introduced where 
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oo — iwt 

X f (w) = _J_ J n(t)X(t)e dt 
2n -oo 

is the Fourier transform of the data as seen through a data 
window n(t). This data window is a real valued function that 
is zero for t<0 and t>T, so that unavailable data are not 
required. And W s is a correction factor, due to the presence 
of the window. 

The estimate S x (w) for a fixed w is a random variable 
with mean 

A 00 00 — iw(tj— 1 2 ) 

E { S x (w) } = W s f dti f dt 2 n (t x ) n (t 2 ) R*(t 2 -tj)e 
47 r 2 ” -oo —oo 

Furthermore, 

^ oo oo -lwT 

E{S X («) } = 1 J W s _ f n(t+r/2)n(t-r/2)dt R x (r)e dr 

27 r —oo 2ir -oo 

with t 1 =t-r/2 and t 2 =t+r/2. Here 

00 

u(r) = W s _ J n (t+r/2) n (t-r/2) dt 
2 tt -oo 

is a lag window u(r) satisfying the following conditions: 

1. u(0) = 1, for preserving power 

A 

2. u(r) = u(-r) which makes S x (w) real 

3 . u (r ) =0 for | t | <T 
For the first condition, 

oo 

u(0) = W s _ J* n z (t) dt=l 

2 7T -oo 


yields 
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W s = 2 

00 

f n 2 (t)dt . 

-00 

Therefore, the estimate becomes 

S x (co) = 2 tt | X f — 

J n 2 (t) dt 
—00 

The second condition, u(r) being even, is obviously satisfied 
while the third condition is also satisfied since u(t) is the 
convolution of two data windows that are only nonzero in the 
range (0,T) . 

The estimate developed above are equivalent in 

expectation to a class of estimates developed by Blackman and 
Tukey in which 

^ oo A -iu >r 

S x (o>) - _1_ Jn(T)-R x (T)e dr 
2n -oo 

/\ 

where R x (r) is an estimate of the autocorrelation function of 
the process and u(r) is a lag window as described above. The 

expectation of this estimate is 

A oo — i o)T 

E (S x (w) ) = _1_ Ju(r)R x (r)e dr. 

27 T -00 

This mean spectral estimate can be shown to equal the 

convolution of the actual spectral density with a "spectral 
window" which is merely the Fourier transform of the lag 
window. Since the autocorrelation is an even function, it can 


be seen that 
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E (S x (w) ) = _±_ Ju(r) ^(r ) coswrdr 
2n 0 

00 

= S x (w) - jR x (r)cos«rdr 

27 r T 

for the case of a "boxcar" lag window, u(r)=l. Therefore, the 
estimate for the power spectral density is biased, but it 
becomes unbiased as T -* In a similar fashion to the first 

technique, we suggest the following technique for estimating 
the two dimensional power spectral density of a nonstationary 


process. The power spectral density for a nonstationary 
process can be written as 

oo oo -i (a; 1 t 1 -a) 2 t 2 ) 

S x («i,u 2 ) - 1 f dt i f dt 2 Rx(t lf t 2 ) e 

air 2 —oo —oo 


= E{X(a> 1 )X*(w 2 ) } , 


where 


00 

X(w) = _1_ J X(t) e~ iut dt . 
2tt -oo 


Here S x (w 1 ,« 2 ) can not be simplified to S x (w) , since the 
autocorrelation R x (t 1 ,t z ) depends on both variables t x and t 2 » 


Therefore , 

from a 

sample 

function of length 

T of a 

nonstationary 

process 

X(t) , 

a class of power 

spectral 


estimates 

S x (wi/W 2 )=W 5 XfC^x) X £ * (u 2 ) 
can be introduced where 

00 -iwt 

X £ («) = _J s n (t) X (t) e dt 

27 T —oo 
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is the Fourier transform of the data as seen through a window 
function n(t) as in the stationary case. 

The mean of the estimate is 

A 00 00 “1 ( CJ^t^ W 2 t 2 ) 

E(S x ( Wl;W2 ) }= W /dt : Jdt 2 rKtJnft^R^t^tJe 

4 7T -CO -00 

The above analysis for the estimate being biased on w x =w 2 in 
the stationary case can be seen to hold when w x is not equal 
to w 2 ■ 

It can be shown [2], that the spectral estimate of the 
finite Fourier transform approach, in the stationary case, is 
essentially a chi-square random variable with two degrees of 
freedom. Therefore, variability of the estimate can be 
reduced by breaking the estimate into N B blocks of 

v\ 

length T B such that N B T B =T . A spectral estimate S xj (w) for 
j=l,2,...,N B will then be taken over each block. If the 
blocks are assumed to be independent, the average of the block 
estimates 


N b a 

S x (c) = S S xj ( w ) 

N b j=l 

is essentially a chi-square random variable with K=2N B degrees 
of freedom. It can be shown with a graph of the variation of 
a chi-square random variable (Figure 2) where a(k) is the left 
bound and b(k) is the right bound that 

P{S x ( w )/a(k) > S x (w) > — S x (w)/b(k) } = 0.80 . 

Therefore, 80 percent of the time the actual spectral density 
will lie between l/a(k) and l/b(k) times the spectral 
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estimate. 

The variability of such estimates is intimately linked 
to their resolution. Full resolution refers to two sinusoids 
of the same amplitude being completely distinguishable when 
viewed through the spectral window function in the frequency 
domain. For the finite Fourier transform technique, full 
resolution requires that frequencies be roughly separated by 

Au> = 27T/T or Af = 1/T . 

Now by breaking the data into N B blocks, the effective length 
has changed from T to T B . Thus, the bandwidth of the estimate 
decreases to 


Af = 1 /T b 

Since, K=2 N b =2T/T b 

K=2AfT 


which shows the tradeoff between variability and frequency 
resolution. 

The estimate in the stationary case 

S x (w) = W s X £ *(w)X f ( W ) 

is essentially a chi-square random variable. Although it is 
not a chi-square random variable, the estimate for the PC 
process 

/■v 

( w i , u>2 ) = W s x £ (fc>i)X £ (w 2 ) 

similarly reduces in variability, when the process is blocked 
in integer multiples n of the period. That is, 

T b = np 

where p is the period of the PC process. 


CHAPTER IV 


SPECTRAL ANALYSIS OF HELICOPTER NOISE 


We will first examine the correlation of the acoustic 
pressure time history X n from an isolated helicopter rotor. 
X n represents a discrete process which reflects the fact that 
the time history is sampled at discrete intervals. It is this 
sampled data which is used to produce a spectral estimate. 
This time history is doubly periodically correlated (DPC) . 
That is 

R(m,n) = EX ra X n * 

is periodic in both m and n. Therefore, its correlation 
function being periodic in m can be written as 

T-l 

(*) R(m, n) = £ R k (n) e 2 " i,nk/T , for all m 

k=0 

And since R(m,n) is also periodic in n, we can write 

R(m,n+T) = R(m,n). 


Therefore, 

T-l T-l 

£ R k (n) e 2 * irak/T = £ R k (n+T)e 2 * i “ k/T , for all m,n 
k=0 k=0 


For a fixed n we have 


T-l 

£ [R k (n+T) - R k (n) ] e 2 " lmk/T = 0 , for all m 
k=0 

which implies that 
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20 


R k (n+T) = R k (n) , for all k 

And since n was arbitrary, R k (n) is periodic for each k. 
Therefore, we can write 


R k (n) 


T— 1 

= Z a>J e-«' T 

j = o 


Substituting into (*) we get 


Hence 


R(m, n) 


T-l T-l 

- 2 [ 2 a ki e 2 " inilT ]e 2 " imk/T . 
k=0 j=0 


T-l T-l 

"D / yy\ \ _ v* V a 2?ri (tnk+n j ) /T 

R(m,n) = 2 2 a kj e 

k=0 j=0 

This means that the spectrum of a DPC process is supported on 

T 2 points { (27rk/T,27rj/T) : j,k = 0,1,..., T-l} (Figure 3). 

Now we can try to get a representation for the process 

X„. We must first utilize the fact that the matrix a kj is 

positive definite. To do this we must show that 

DAD* = 2 d k a k jdj > 0 
k, j 

for each vector D=(d 0 ,d 1 , . . . e C T . Let (c 0 , c : , . . . , c T . x ) be 

a vector in C T such that 


T-l 

2 C n e 2 * ink/T = d k , for every k=0, 1, . . . ,T-1. 
n=0 


This is possible by Fourier theory. Now we can write 

T-l 


T-l 
2 d k a 
k, j=0 


•kj d j* = 2 


T-l T-l 

(2 c ro e 
k,j=0 m=0 


2nimk/T 


) a k j (2 c n e 
n=0 


2ninj /T 


T-l 
= 2 
m=0 


T-l T-l 
E c m 
n=0 k,j=0 


(2 a kj e 


2ni(mk-nj )/T 


)C„^ 
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T-l T-l 

= S E c m R (m, n) c n * 
m=0 n=0 

T-l T-l 

= E{(S cAHS c„,X„)*} 
m=0 m=0 


= || Z c n X n || 2 > 0 Q.E.D. 

n=0 

Since a k j has been shown to be positive definite, we know there 
is a vector (y 0 ,y 1 , . . . ,yj-j) of Gaussian random variables with 
mean (0,0,...,0) and covariance matrix (a kj ) . So 

a kj = E(y k yj*) for every k,j = 0,1, — ,T-1. 
Therefore, we can write (4.2) as 


R x (m,n) = Z S E(y k yj)e 2 " i(M)/T 
k=0 j=0 

= E{^2 1 e 2 ’ ,imk/T y k ) (Se 2,rinJ/T yj )*} 
k=0 j — 0 

This means that 

T-l 

X n = E e 2 * ink/T y k , for every n=0 , ±1 , ±2 , . . . 
k=0 

Gi v - i . » 2nin/T , _(2)2*in/T , . (T-l)27rin/T 

X n = y 0 + y x e + y 2 e +...+y T _ 1 e 


Therefore, 

X n = X n+T . 

And since X n is periodic, prediction is easy once the period 
of the process is known. 

The object of this thesis has been to develop a technique 
for analyzing helicopter noise in a more exact way. The 
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standard approach is to treat the data as if it were 
stationary and to produce a one dimensional spectrum from 
which the data will be analyzed. In actuality the data 
produced from the isolated tail rotor and isolated main rotor 
is periodically correlated which requires a two-dimensional 
spectrum to completely represent the spectrum of the data. 
Further the data from the main and tail rotor combination is 
basically the sum of two incommensurate periodic components 
which means that this data is neither periodic nor stationary. 
Hence, this data also requires a two dimensional spectrum. 

In tryinq to develop a technique to analyze helicopter 
noise the necessity for the spectrum to be two-dimensional is 
utilized. As discussed in chapter 3, an estimate was 
developed for the power spectrum of a nonstationary process. 
And it was this estimate that was used to produce two 
dimensional spectral estimates from the data. The technique 
consists of estimating the spectrum over a region of values 
for wj and o> 2 and then block averaging the estimate to reduce 
variability. In the case of the isolated main and tail rotors 
blocking was done as a multiple of the period of process. 
This was done because of the characteristics of a periodic 
process. A periodic process has an amplitude spectrum as 
discussed above with support only on lattice points. And to 
insure that these points are actually viewed, a length of data 
which is a multiple of the period of the process must be used 
when the data is Fourier transformed. 
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The data used to produce the spectral estimates in this 
thesis was taken from a wind tunnel test conducted at 
NASA/Langley. The Sikorsky Aircraft's Basic Model Test Rig 
was used in the experiment details of which can be found in 
[6]. The data was taken from several locations around the 
helicopter model (see Figure 4) . By analyzing data from 
several different positions, the noise pattern produced could 
be used to determine how the noise is radiated and in which 
direction most of the noise travels. 

In producing the two dimensional spectrum for the 
isolated tail rotor, program 1 (see Appendix A) was created 
to produce values of the estimate. The tail rotor had a blade 
passage frequency of around 450 Hz which was determined by 
graphed data (see Figure 5, lines for the shaft frequency are 
also present) produced by computer program 2 (see Appendix B) . 
The motivation for this program was discussed in chapter II. 
The blade passage frequency was also verified analytically 
using the given test conditions. The passing blades produce 
a corresponding periodic acoustic pressure time history with 
a fundamental frequency of 450 Hz. It is this time history 
from which the sampled data was taken. After the data was 
input to program 1 (Appendix A) , a set of values were given 
in the output for the spectral estimate at frequencies wj. and 
w 2 going from zero to a chosen upper frequency (for us 4800 
Hz) . We viewed the data and decided that due to the 
background noise at low frequencies the scaling of the graph 
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did not allow enough of the power at the fundamental frequency 
and at harmonics of the tail rotor period to be shown. 
Therefore, we decided to filter the data using a Chebyshev 
digital filter. All frequencies below 250 Hz were removed, 
and a graph of the output data was viewed once again. This 
time more detail of the power at the fundamental frequencies 
and its harmonics could be seen (see Figures 6, 7, 8, 9). 

The graphs show most of the power for the tail rotor at 
the fundamental frequency (450 Hz) and the first and second 
harmonics (900 and 1350 Hz respectively) . There are also 
components off the main diagonal which represent the 
correlation of the amplitude of sinusoids at frequencies 
and u> 2 where * u> 2 » 

As for the data sampled from the isolated main rotor time 
history, this data is periodic with a frequency of 95 Hz. The 
first graphs of this data showed high power levels at around 
95 Hz. Since the background noise was also high around 95 Hz, 
we decided to filter frequencies below 250 Hz out of this data 
also. This resulted in the second harmonic of the main rotor 
noise (285 Hz) being dominant. At microphone 3 there is 
noticeable power at the third and fourth harmonics also. At 
microphone 5 the background tonal noise between 800 and 1300 
Hz as reported in [6] is apparent. See figures 10, 11, 12 and 
13 for the spectral estimates produced from data taken by 
microphones 2 through 5. 

The spectrum of the combined noise was of interest due 
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to the fact that if noise for the isolated main and tail 
rotors is independent, the theory shows that the spectrum for 
the combined noise is just the sum of the spectra of the 
isolated main and tail rotors. After filtering the low 
frequencies as above, it could be seen from figures 14, 15 ,16 
and 17 that to some degree the combined noise had a spectrum 
which was the sum of the main and tail rotor spectra at 
corresponding microphones. However, the question was raised 
as to the independence of the noise from the main and tail 
rotors when both rotors are operating due to the possibility 
of the main rotor wake being swept back into the tail rotor. 

We now took another approach. We took the data from the 
main rotor and added it directly to the data from the tail 
rotor. Then we produced the two dimensional spectrum (see 
Figures 18, 19, 2 0 and 21) for this added noise. As the 
theory shows, these spectral estimates are the sum of the 
spectra of the isolated main and tail rotors at corresponding 
microphones. After comparing spectrums of the added noise to 
those of the combined noise, we concluded that there is some 
degree of dependence between the noise from the main rotor and 
tail rotor when the combined noise is being produced. 

It should be noted again that the combined noise is not 
periodically correlated. It is the sum of two periodically 
correlated processes with incommensurate periods. According 
to theory, if the two periodically correlated processes are 
independent, the resulting spectrum should equal the sum of 
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the spectra of each process. Since the combined noise can 
only be blocked with respect to one period, we chose to block 
the data with respect to the period of the tail rotor noise. 
This will naturally result in the estimated spectrum being 
skewed to some degree. However, the same amount of skewing 
occurs in the estimated spectrum for the added noise which is 
also blocked with respect to the tail rotor noise. Therefore, 
the estimated spectrum of the combined and added noise can be 
viewed equally with respect to theory. 

A. Assumed Stationaritv 

The standard method for handling data taken from 
helicopter noise is to treat the data as if it were 
stationary. This results in a one dimensional spectrum along 
the u>i=io 2 diagonal. Therefore, the spectrum excludes all 
information concerning correlation of the sinusoidal 
amplitudes at frequencies where Although these 
correlations need further investigation, they can not be 
neglected if the spectrum of a nonstationary process is to be 
completely studied. Examples of sound pressure spectra, for 
the tail rotor produced using this method, are given in 
figures 22, 23, 24 and 25. The graphs produced contain 22 
degrees of freedom and a frequency resolution of 49 Hz. The 
noticeable peaks are the fundamental frequencies and 
subsequent harmonics. In the graphs the pressure level drops 
off very rapidly. Notable peaks start to drop off after 3000 
Hz. Due to the frequency range of audible sound (20 hz to 20 
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khz) , the ear suggests high tonal levels above 3000 Hz also. 
However, this is not evident with this technique. 

B. Phase Shifted Data 

A periodically correlated process will become stationary 
by applying a random phase shift which is uniformly 
distributed over the period of the process [7]. This 
stationary process can then be adequately analyzed by the 
standard method discussed above. This new stationary process 
which has been produced is no longer the original process. 
And thus, no longer contains information about the correlation 
of Fourier components at different frequencies, Hardin and 
Miamee [3]. To implement the technique, we take a length of 
data and break it into blocks. Each block is of length 
necessary to obtain a desired frequency resolution plus the 
period of the underlying process. When the program is 
implemented, a random function call chooses a sample index 
uniformly distributed over the period of the process. 
Starting with this value, the program Fourier transforms 
enough samples to produce the desired frequency resolution. 
This is done for each block of data (see program 3 Appendix 
C) . This technique results in wastage of sample values from 
the start of the record up to the chosen index, but does 
implement the random phase shift while maintaining the desired 
number of samples per block. 

The sound pressure spectra produced from this shifted 
data cam be compared to that of the unshifted data one 



28 


dimensional spectra. After viewing the graphs in figures 
26,27, 28 and 29, it is apparent that there is significant 

power at frequencies above 3000 Hz. The impression which we 
get when overlaying graphs of the shifted and unshifted data 
is that by not shifting the data, we get a smoothed spectrum 
as if a moving average had been applied to the spectrum. 
There are additional peaks in the spectrum which appear to 
correspond to the shaft frequency. Since the shaft frequency 
is one fourth of the blade passage frequency, there could be 
three peaks between the harmonics of the spectrum for the tail 
rotor noise. This is evident in [6] in which a frequency 
resolution of 12 Hz is used. This resolution is fine enough 
to show the shaft harmonics. However, our graphs have a 
frequency resolution of 50 Hz which is not fine enough to 
completely show the shaft harmonics. This technique of 
shifting the data appears to be useful in resolving harmonics 
of the transformed data, and also rests on a firm theoretical 
foundation. 

The output resulting from this method was also used in 
order to produce a two dimensional spectral estimate. The 
graphs for the data taken at microphones 2 through 5 are given 
in figures 30, 31, 32 and 33. The graphs contain off diagonal 
values which are not to be expected for stationary data. We 
suspect that this is due to our limited data length. The time 
history could only be broken into eleven blocks in which the 
data is shifted. This small amount of data does not suffice 
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for a uniform representation over the period of the process, 
which is required in order to produce a stationary process 
from a PC process. 

CONCLUSION AND RECOMMENDATIONS 

In conclusion we feel that we have successfully developed 
a method for producing a two dimensional spectral estimate. 
This estimate is known to contain more information than the 
one dimensional estimate which results from considering 
nonstationary data as being stationary. We have shown that 
applying a random shift to a periodically correlated process 
results in a more useful spectrum for viewing the higher 
harmonics when the data is analyzed as stationary. The 
harmonics resulting from the standard method discussed above 
appear to drop off to quickly. 

For further study, we recommend that values off the main 
diagonal of the spectral estimate be investigated as to how 
their presence can be used to characterize the process itself. 
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Figure 6. Tail Rotor Spectrum Microphone 2 
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Figure 8. Tail Rotor Spectrum Microphone 4 
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Figure 9 . Tail Rotor Spectrum Microphone 5 
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Figure 10. Main Rotor Spectrum Microphone 2 
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Figure . Combined Spectrum Microphone 2 
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Figure 19. Added Spectrum Microphone 3 
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Figure 20. Added Spectrum Microphone 4 
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Figure 22. Tail Rotor Sound Pressure Level 
Microphone 2 
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Figure 23. Tail Rotor Sound Pressure Level 
Microphone 3 
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Figure 24. Tail Rotor Sound Pressure Level 
Microphone 4 
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Figure 25. Tail Rotor Sound Pressure Level 
Microphone 5 
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Figure 26. Tail Rotor Sound Pressure Level 
with Random Shift Microphone 2 
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Figure 27. Tail Rotor Sound Pressure Level 
with Random Shift Microphone3 
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Figure 28. Tail Rotor Sound Pressure Level 
with Random Shift Microphone 4 
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Figure 29. Tail Rotor Sound Preessure Level 
with Random Shift Microphone 5 
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Figure 30. Tail Rotor Spectrum with Random Shift 
Microphone 2 
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Figure 31. Tail Rotor Spectrum with Random Shift 
Microphone 3 
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Figure 32. Tail Rotor Spectrum with Random Shift 
Microphone 4 
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Figure 33. Tail Rotor Spectrum with Random Shift 
Microphone 5 
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C234567 

INTEGER N,IWK(6150) ,S,Q 
REAL A(840) ,WK(6150) 

DIMENSION P(IOOOO) 

COMPLEX XWT(0 : 500 ,11) , SXW, SUM, SUMl 
COMPLEX X(1010) 

B-840 

S-840 

DO 6 N“1 , 10000 
READ (1 , *) P(n) 

6 CONTINUE 

Q--839 

DO 40 J-1,11 

Q-Q+340 

K-0 

DO 8 I-Q.Q+839 

K-Ktl 

A(K)-P(I) 

8 CONTINUE 

CALL FFTRC (A, S ,X, IWK.WK) 

ND2-S/2 
DO 20 1-2 , ND2 

X(S+2 - 1)“C0NJC(X( I ) ) 

20 CONTINUE 

DO 50 1-1 , li/2 + 1 
XWT ( I - 1 , J ) -X ( I ) 

50 CONTINUE 
40 CONTINUE 

DO 90 Q-0,100 
DO 101 M-0,100 
SUM“(0. 0,0.0) 

DO 70 J-1,11 

SUM- . 0011108*(CONJC(XUT(M, J) )*XWT(Q , J) )+SUM 
70 CONTINUE 

SUM1-SUM/11 . 0 
sxw-CABS (SUMl ) 

WRITE ( 6 , *) M , Q , SXU 
101 CONTINUE 
90 CONTINUE 
END 


Appendix A. 2-D Spectrum 
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INTECER N , IWK(10152 ) , S 
REAL SUM1 , SUM2 , SUM3 
REAL A(3334) ,WK(10152) 

DIMENSION U(10000) 

COMPLEX X(0: 1668) , SUM 
S-3334 

DO G N-1,3334 
READ( 1 , *)U(N) 

6 CONTINUE 

DO 8 1-1,3334 

A ( I ) **U ( I ) 

8 CONTINUE 

CALL FFTRC (A , S , X , IWK , WK) 

ND2-S/2 
DO 20 I-1.ND2 

X(S+2 - 1 ) -CONJC (X( I ) ) 

20 CONTINUE 

DO 7 1-1, S/2+1 
X(I-1)-X(I) 

7 CONTINUE 

DO 200 Q-0, 100,1 
DO 100 P-0,100,1 
SUM-(0. 0,0.0) 

DO 10 M-0,7 

SUM-SUM + X( P+M)*CONJ C (X(Q+M) ) 

10 CONTINUE 

SUMi-( (REAL(SUM) )**2)+( (AIMAC(SUM) )**2) 

SUM2-0.0 
DO 21 M-0,7 

SUM2-SUM2 +( (REAL(X(P+M) ) )**2) +( (AIMAG(X(P+M) ) )**2) 

21 CONTINUE 
SUM3-0 . 0 

DO 30 M-0 , 7 

SUM3-SUM3 + ( (REAL(X(Q+M) ) )**2) +( (AlMAG(X(Q+M) ) )**2) 
30 CONTINUE 

F-SUM1/ ( SUM2*SUM3 ) 

C-.3 3 

R-P*(20000. 0/1667 . 0) 

U-Q*( 20000 . 0/1667 . 0) 
if (F.ge.C)then 
URITE(6 ,*)R,W 
ENDIF 

100 CONTINUE 
200 CONTINUE 
END 


Appendix B. Spectral Coherence 


66 



C234567 

INTECER N,IWK(6150) ,S,Q 
REAL A(820) ,WK(6150) 

DIMENSION P(10000) 

COMPLEX XUT(0 : 500 , 11) , SX, SUM, SUM1 
COMPLEX X(1010) 

B-820 

S-820 

DO 6 N-l, 10000 
READ(l,*)P(n) 

C CONTINUE 
Q--907 

DO 40 J-1,11 
Q-Q+908 
K-RAN(Q)*89 
N-0 

. DO 8 I-Q.Q+907 
N-N+l 

IF(N . LE . 820)THEN 

K-K+l 

A(N)“P(K) 

ENDIF 

8 CONTINUE 

CALL FFTRC (A , S , X , IWK, WK) 

ND2-S/2 
DO 20 I“2 , ND2 

X(S+2-l) — CONJC (X( I ) ) 

20 CONTINUE 

DO 50 I-l.B/2+1 
XWT (I-1,J)-X(I) 

50 CONTINUE 
40 CONTINUE 

DO 101 M-0,100 
SUM-(0. 0,0.0) 

DO 70 J-1,11 

SUM- . 0011108*(CONJC(XWT(M, J) )*XWT(M, J) )+SUM 
70 CONTINUE 
SX-SUM/11.0 
SXW-REAL(SX) 

G-SXW 

SPL-10*LOC10 (C) + 74 
R-M* (20000/4 10) 

WRITE (6 , *)R , SPL 
101 CONTINUE 

run 


Appendix C. l-D Spectrum 
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ABSTRACT 


A TWO DIMENSIONAL POWER SPECTRAL ESTIMATE 
FOR SOME NONSTATIONARY PROCESSES 

Student: Gregory L. Smith Department: Mathematics 

Date of Birth: July 6, 1961 Advisor: Dr. Abol Ghassem Miamee 

A two dimensional spectral estimate for a nonstationary 
process is developed. It addresses the need for modelling 
phenomena which exhibit nonstationary behavior. Currently 
these phenomena are usually modelled as stationary processes. 
The spectrum of a nonstationary process is two dimensional 
while the spectrum of stationary processes is one dimensional 
Therefore, these phenomena are not as completely represented 
when modelled as stationary processes. The usual method of 
analyzing nonstationary processes as if they were stationary 
is compared to the two dimensional estimate which has been 
developed. In addition, a random phase shift which as the 
theory shows should produce a stationary process, is 
introduced to our nonstationary process. This stationary 
process can then be analyzed in the usual way. The results 
are compared with the methods described above. 



